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During  the  past  year  we  have  initiated  a  new  research  program  to  experimentally  and 
theoretically  investigate  ultrafast  processes  and  carrier  dynamics  in  compound  semiconduc¬ 
tors.  The  objective  of  our  program  is  to  apply  state  of  the  art  femtosecond  measurement 
techniques  including  high  resolution  pump-probe  absorption  spectroscopy  and  time  division 
interferometry  measurements  of  nonlinear  index  as  well  as  advanced  theoretical  techniques 
including  ensemble  Monte  Carlo  calculations  and  analytic  solutions  of  rate  equation  mod¬ 
els  to  study  transient  processes  in  semiconductors.  The  combination  of  experimental  and 
theoretical  approaches  can  be  used  to  provide  fundamental  information  about  the  physics 
of  excited  carriers  in  semiconductors  and  how  they  impact  on  electronic  and  optoelectronic 
device  performance. 

Experimental  efforts  at  MIT  have  centered  on  the  development  of  new  techniques  for 
measurement  of  ultrafast  processes  in  semiconductors  and  their  application  to  study  tran¬ 
sient  processes  in  AlGaAs.  Studies  focus  on  both  femtosecond  measurements  of  nonlinear 
index  as  well  as  absorption.  Investigations  of  nonlinear  index  are  relevant  for  applications 
in  all  optical  switching  as  well  as  high  speed  modulation  using  semiconductor  devices. 
During  the  past  contract  year,  we  have  developed  a  new  femtosecond  time  division  inter¬ 
ferometry  technique  for  characterizing  nonlinear  index  changes  from  different  components 
of  the  tensor.  This  new  technique  uses  pump  and  probe  measurements  to  investigate 
the  amplitude  and  femtosecond  dynamics  of  the  nonlinear  phase  shift  which  is  induced 
by  a  pump  pulse  on  a  probe  pulse.  Because  nonlinear  index  effects  are  associated  with 
below  bandgap  excitation,  they  are  extremely  small,  and  sensitive  measurements  using 


92-12577 

iiiiiiiniiiiii 


92  5  II  049 


1 


waveguide  geometries  are  required.  We  have  developed  a  new  time  division  multiplexed 
interferometer  scheme  which  permits  the  measurement  of  nonlinear  phase  shift  without 
parasitic  effects  from  thermal  index  changes  or  acoustic  vibrations.  This  new  technique 
permits  the  characterization  of  the  parallel  and  perpendicular  components  of  the  nonlinear 
index.  Preliminary  studies  have  been  performed  using  optical  fibers  and  semiconductor 
waveguides.  Future  investigations  will  apply  this  technique  to  characterize  nonlinear  index 
in  compound  semiconductors  as  well  a s  quantum  confined  structures.  Future  studies  will 
include  investigations  in  the  strained  layer  InGaAs/GaAs  systems. 

Other  ongoing  experimental  work  at  MIT  focuses  on  studies  of  femtosecond  carrier 
dynamics  in  AlGaAs.  We  are  currently  performing  pump  probe  measurements  of  carrier 
dynamics  in  AlGaAs  using  tunable  40  fs  pulses  from  a  high  repetition  rate  femtosecond 
amplifier.  Pump  probe  absorption  saturation  measurements  are  an  indicator  of  the  scat¬ 
tering  carriers  from  their  initial  optically  excited  states.  By  varying  the  wavelength  of 
the  excitation  pulses,  the  excess  energy  of  the  carrier  distributions  may  be  continuously 
varied.  The  objective  of  these  studies  is  to  investigate  intervalley  scattering  by  varying 
the  carrier  energies  above  and  below  the  allowed  energy  for  scattering  from  the  T  to  the 
L  valleys. 

Finally  we  are  beginning  preliminary  work  on  three  pulse  experiments  to  investigate 
the  effect  of  a  cold  carrier  distribution  on  the  relaxation  processes  of  hot  carriers.  In 
these  measurements  an  initial  electron  and  hole  distribution  is  prepared  by  exciting  the 
sample  with  a  variable  intensity  ultrashort  pulse.  After  a  suitable  delay  to  permit  the 
excited  characters  to  thermalize,  a  standard  pump  probe  measurement  is  performed.  Pre¬ 
liminary  measurements  have  been  performed  using  2.0  eV  pulses  of  40  fs  duration  from 
a  colliding  pulse  modelocked  laser  source.  Future  experiments  using  tunable  pump  and 
probe  measurements  will  permit  the  investigation  of  both  the  scattering  of  excited  carriers 
with  a  cold  carrier  background  as  well  as  the  effect  of  high  energy  carrier  excitation  on 
a  quasi-equilibrated  electron  in  whole  distribution.  These  experiments  have  direct  bear¬ 
ing  on  gain  and  absorption  dynamics  in  laser  diodes  where  perturbations  of  quasithermal 
carrier  processes  produce  changes  in  gain  and  index  in  high  speed  modulation. 

Theoretical  work  at  the  University  of  Florida  has  focused  on  developing  new  approaches 
for  studying  ultrafast  carrier  dynamics  of  carriers  photoexcited  in  GaAs  and  AlGaAs.  In¬ 
vestigations  have  concentrated  on  Monte  Carlo  calculations  including  both  the  effects  of 
electron  and  hole  distributions  and  carrier-carrier  scattering  as  well  as  analytic  rate  equa¬ 
tion  models  of  pump/continuum  probe  and  tunable  pump  probe  spectroscopy.  Because 
the  dynamics  of  excited  carriers  in  semiconductors  is  extremely  complex,  detailed  theoret¬ 
ical  studies  such  as  ensemble  Monte  Carlo  or  rate  equation  models  are  required  to  extract 
fundamental  information  on  carrier  dynamics  from  experimental  measurements.  During 
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the  past  year,  we  have  performed  a  number  of  theoretical  investigations  which  are  aimed 
at  unifying  experiment  and  theory. 

We  have  developed  ensemble  Monte  Carlo  calculational  techniques  and  correlated  pre¬ 
dicted  measurements  of  pump  and  probe  absorption  saturation  spectroscopy  in  GaAs  and 
AlGaAs.  These  studies  resulted  in  the  identification  and  calculation  of  Coulomb  enhance¬ 
ment  effects  from  excitons  near  the  band  edge. 

Other  Monte  Carlo  studies  have  shown  the  importance  of  carrier  diffusion  effects  in 
interpreting  the  role  of  carrier  dynamics  in  both  semiconductors.  Theoretical  investiga¬ 
tions  demonstrated  that  carrier  diffusion  can  significantly  effect  the  results  of  pump  probe 
absorption  saturation  measurements  of  carrier-carrier  scattering  on  the  time  scale  of  10- 
100  fs.  These  investigations  demonstrate  that  differences  in  sample  thickness,  even  on  the 
scale  of  0.2-0. 5  /mi,  affect  the  interpretation  of  experimental  results  by  altering  the  carrier 
density  profiles  in  the  semiconductors. 

Finally,  we  have  developed  an  analytic  solution  of  a  rate  equation  model  to  describe 
intervalley  scattering.  The  rate  equation  approach  provides  a  simpler  and  computation¬ 
ally  less  complex  model  for  interpreting  femtosecond  absorption  saturation  spectroscopic 
results.  This  model  has  a  direct  bearing  on  a  number  of  previous  experiments  which  at¬ 
tempt  to  measure  intervalley  scattering  rates.  Solutions  of  a  rate  equation  approach  show 
that  the  return  time  of  electrons  from  the  satellite  L  and  X  valleys  is  limited  not  by  the 
intervalley  scattering  rate  but  is  instead  limited  by  inelastic  scattering  within  the  T  valley 
(chiefly  through  polar  optical  phonon,  POP  emisssion).  These  results  suggest  that  the 
accurate  measurement  of  intervalley  scattering  must  depend  on  the  measurement  of  the 
initial  scattering  rate  from  T  to  L  and  not  the  return  rate  by  which  carriers  return  to  the 
T  valley. 

Future  theoretical  studies  will  extend  these  results  and  develop  a  closer  correspondence 
between  ensemble  Monte  Carlo  and  rate  equation  models  and  femtosecond  spectroscopic 
studies.  The  objective  of  these  investigations  is  to  extract  fundamental  information  on 
carrier  scattering  events  such  as  intervalley  scattering  and  role  of  carrier-carrier  scattering 
including  electron  hole  distributions.  Additional  studies  will  extend  these  results  to  inves¬ 
tigate  coherent  and  non-Markovian  processes,  i.e.,  quantum  corrections  to  the  standard 
Boltzmann  theory  of  carrier  relaxation,  at  very  short  time  scales  of  less  than  15  ps. 
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Femtosecond  time  division  interferometry  technique  for  measuring 
the  tensor  components  of  *(3) 

C.  de  C.  Chamon,  C.  K.  Sun,a)  H.  A.  Haus,  and  J.  G.  Fujimoto 

Department  of  Electrical  Engineering  and  Computer  Science.  Research  Laboratory  of  Electronics, 
Massachusetts  Institute  of  Technology.  Cambridge,  Massachusetts  02139 

(Received  3  September  1991;  accepted  for  publication  16  November  1991) 

We  describe  a  new  femtosecond  time  division  interferometry  technique  for  characterizing 
nonlinear  index  changes  from  different  components  of  the  ^,3)  tensor.  Pump  probe 
interferometric  measurements  of  nonlinear  phase  are  performed  using  a  time  division 
multiplexed  reference  pulse.  The  /t2||  and  tin  components  of  the  nonlinear  index  are  measured 
by  permuting  the  polarizations  of  the  pulses  and  using  a  novel  differential  detection  and 
modulation  scheme.  Measurements  are  demonstrated  in  an  optical  fiber. 


Nonlinear  index  changes  have  been  investigated  using 
a  number  of  experimental  techniques  including  four-wave 
mixing,1  nonlinear  waveguide  couplers,2'3  nonlinear 
Fabry-Perot,4  fringe  shift  interferometry,5  Mach-Zehnder 
interferometry,6  and  beam  propagation  distortion.7  In  gen¬ 
eral,  direct  measurements  of  index  nonlinearities  are  com¬ 
plicated  because  the  measurement  techniques  are  sensitive 
to  thermal  or  acoustical  parasitic  signals,  have  low  sensi¬ 
tivity,  or  require  deconvolution.8 

Recently,  we  have  developed  a  new  technique,  femto¬ 
second  time  division  interferometry  (TDI),  that  permits 
high-sensitivity,  direct  measurements  of  the  nonlinear  in¬ 
dex  n2  in  waveguide  devices.4,10  TDI  was  previously  lim¬ 
ited  to  measurements  of  n2i,  i.e.,  index  changes  induced  on 
a  probe  pulse  by  an  orthogonally  polarized  pump.  How¬ 
ever,  new  materials  such  as  organic  polymers,  multiple 
quantum  wells  and  quantum  wires  are  anisotropic  and  a 
complete  characterization  of  n2  requires  the  measurement 
Of 

Time-resolved  measurements  of  n2||  are  complicated  by 
the  fact  that  both  the  pump  and  probe  pulses  have  the  same 
frequency  and  polarization.  In  TDI  measurements,  this  de¬ 
generacy  prevents  the  interference  signal  between  the 
probe  and  reference  from  being  isolated  from  the  interfer¬ 
ence  signals  between  the  pump  and  probe  or  the  pump  and 
reference.  In  this  letter  we  report  a  new  extension  of  the 
time  division  interferometry  technique  that  permits  mea¬ 
surements  of  «2||.  The  contributions  to  n2  from  different 
components  of  ^<3)  can  be  characterized  by  measuring  dif¬ 
ferent  permutations  of  pump  and  probe  polarizations. 

The  laser  source  used  for  our  experiments  was  a  mode- 
locked  Nd:YAG  which  was  pulse  compressed  and  fre¬ 
quency  doubled  to  synchronously  pump  a  dye  laser.  Using 
the  dye  Styrl  9,  the  system  generates  subpicosecond  pulses 
tunable  from  795  to  875  nm.  The  experimental  configura¬ 
tion  for  time  division  interferometry  has  been  described 
previously  .9 10  A  portion  of  the  laser  output  is  directed  to  a 
micrometer  stage  delay  line  to  function  as  the  pump.  The 
pump  beam  is  chopped  and  the  phase  delay  of  the  pulse  is 
modulated  using  a  PZT  attached  to  one  of  the  delay  mir¬ 
rors.  The  probe  and  reference  are  obtained  using  a  Mich- 
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elson-type  delay  line  with  waveplates  and  a  polarizing 
beam  splitter.  This  system  functions  as  a  polarization  sen¬ 
sitive  delay  line  and  generates  orthogonally  polarized 
probe  and  reference  pulses.  Half-wave  plates  placed  in  the 
probe,  reference,  and  pump  beams  are  used  to  set  different 
combinations  of  polarizations  for  the  three  pulses.  The 
pulses  are  then  coupled  into  the  waveguide  structure  to  be 
characterized.  After  the  waveguide,  the  signals  are  de¬ 
tected  using  a  polarizer  and  dual  detector.  Finally,  the 
phase  bias  operating  point  of  the  interferometer  is  actively 
stabilized  using  a  PZT  on  the  polarization-sensitive  delay 
line.  The  feedback  error  signal  is  obtained  using  a  boxcar 
integrator  which  detects  the  interferometer  output  signal 
when  the  pump  pulse  is  blocked. 

Figure  1  shows  the  relative  positioning  of  the  pulses  for 
the  measuring  n2|i  a°d  nu-  The  field  amplitudes  of  the 
pump,  probe,  and  reference  pulses  are  denoted  Epu,  Epr, 
and  EKt-  The  pump  pulse  induces  a  nonlinear  phase  shift 
on  the  probe  which  is  then  measured  by  interfering  the 
probe  pulse  with  a  time  division  multiplexed  reference 
pulse.  The  reference  pulse  is  then  time  delayed  using  a 
polarization-sensitive  delay  line  so  that  it  interferes  with 
the  probe  pulse.  For  the  measurement  of  n2||  the  pump 
pulse  does  not  have  a  different  polarization  and  thus  is  not 
separated  from  the  probe  and  reference  pulses.  The  phase 
of  the  interferometer  can  be  biased  by  adjusting  the  polar¬ 
ization-sensitive  delay  line  with  a  PZT.  The  nonlinear  in¬ 
dex  is  characterized  by  measuring  the  phase  shift  of  the 
probe  relative  to  the  reference  using  a  polarizer  and  differ¬ 
ential  detector  geometry.  Finally,  the  dynamics  of  the  non¬ 
linear  index  may  be  measured  by  varying  the  delay  be¬ 
tween  the  pump  and  probe  pulses. 

The  details  of  the  measurement  technique  for  n2||  and 
n2l  can  be  described  by  performing  a  simple  analysis  of  the 
signals  which  are  detected.  Let  us  first  consider  the  case 
where  the  pump  and  probe  are  perpendicularly  polarized. 
The  signals  projected  along  the  two  polarization  axes  of  the 
detectors  are 

51.2~|£p^"l±£rri^|2+  If^l2,  (1) 

where  <t>„\  is  the  nonlinear  phase  shift  produced  by  the 
pump  on  the  probe,  and  is  the  bias  phase  of  the  inter¬ 
ferometer.  In  order  to  simplify  the  equations,  we  express 
the  field  amplitudes  as  real  variables  with  relative  phases 
expressed  explicitly.  The  square  of  (he  field  amplitude  is 
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FIG.  1  Relative  positioning  of  the  pump,  probe,  and  reference  pulses 
before  and  after  the  waveguide  sample  for  measurements  of  n.  and  » 
The  polarization  sensitive  delay  line,  D.  separates  pump  and  probe  in  the 
case  of  the  perpendicular  measurement,  but  not  in  the  parallel  case 
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normalized  to  intensity.  Performing  differential  detection 
cancels  the  background  intensity  from  probe  and  reference 
fields  and  yields  the  interference  term: 

S~4£pr£rt.f  cos(4,  -  4,as).  (2) 

For  measurements  of  nonlinear  phase,  the  interferometer  is 
biased  to  <fShias  =  tt/2,  and  the  phase  changes  «Snl  are  small 
so  the  signal  is  given  by 

S  ~  4 £pr£rcf  sin  4,  ~  4£pr£rr(4i'  ( 3 ) 

Since  the  pump  beam  is  chopped  and  4i  is  modulated, 
lock-in  amplification  detects  only  the  terms  that  con'ain 
£pu.  The  contribution  of  these  terms  to  the  nonlinear  phase 
shift  is 


1  -f-  cos  2 <t> 


X(,2J2(<u:6),  -  (o,(o ) 
2Tl22l(<‘,:w>  ~ 


where  4>  is  the  phase  between  pump  and  probe  and  nn  is 
defined  as  the  Kerr  coefficient  for  monochromatic  orthog¬ 
onal  polarizations  with  the  coherence  term  suppressed. 
The  term  containing  cos  2 i>  is  the  result  of  cross-phase 
modulation  between  the  pump  and  reference.  This  term  is 
averaged  to  zero  when  the  pump  phase  is  dithered  using  a 
PZT,  or  when  the  pump  time  delay  is  scanned.  The  result¬ 
ing  signal  permits  a  measurement  of  n21. 


S~4EprEnfElu(2n/A)nnLt(r.  (5) 


Next  let  us  consider  the  case  where  the  pump  and  probe 
are  polarized  parallel  to  each  other  and  interfere.  For  this 
case  the  detected  signals  in  the  two  polarizations  are 


S|.2~l(£pr  +  )  e*"  ±  Ert(  e*”«|2,  ( 6 ) 

where  <b  is  the  phase  difference  between  probe  and  pump. 
Differentia]  detection  yields 


FIG.  2  Time  division  interferometer  measurement  performed  with  and 
without  dithering  the  phase  of  the  pump  pulse. 

S ~ 4 £pr£ref cos ( 4,  -  4, *. ) 

+  4£pu£rcr cos ( d>  4  4,  -  4,hs)-  (2 ) 

Biasing  at  4ias  =  f/2  for  small  <t>nl  the  detected  signal  is 

S~4£pr£r(.f  sin  4,  4-  4£pu£rcf  sin(<b  +  4,) 

~ 4 £pr£rc,4,  +  4£pu£rcf( sin  6  +  4i  cos  6).  ( 8 ) 

The  second  term  in  (7)  and  (8)  arises  from  interference 
between  the  pump  and  reference  pulses. 

The  nonlinear  phase  shift  produced  by  the  pump  on 
the  probe  is 

4i  =  ( 2 tr/A ) n2,| ( £2r  +  £pu  +  2£pr£pu  cos  <t)Lrff,  ( 9 ) 

where  n2(|  is  the  Kerr  coefficient  for  a  monochromatic  lin¬ 
early  polarized  wave.  The  third  term  in  (9)  is  the  interfer¬ 
ence  of  the  pump  and  probe  fields. 

In  order  to  suppress  the  phase-sensitive  cross  terms 
and  perform  a  measurement  of  n2 y,  the  delay  of  the  pump 
pulse  is  dithered  by  an  integer  multiple  of  the  wavelength. 
This  averages  the  terms  containing  sin  <b  and  cos  <i>  to  zero 
and  the  term  containing  cos2  <j>  to  1/2.  Since  lock-in  detec¬ 
tion  is  used,  the  chopping  frequency  of  the  pump  must  be 
lower  than  the  frequency  of  the  pump  phase  to  average  the 
phase-sensitive  contributions  in  the  detected  signal.  After 
lock-in  amplification  the  detected  signal  is 

S~  8£pr£rer£^u(2^//l  )nlfLen.  ( 10) 

Note  that  the  detected  signal  for  the  parallel  case  has  ex¬ 
actly  the  same  form  as  in  the  perpendicular  case  except  for 
an  extra  factor  of  2. 


534 


Appl  Phys.  Lett ,  Vot  60,  No  5,  3  February  1992 


Chamon  el  al. 


534 


FIG.  3.  Time  division  interferometer  measurement  of  the  nonlinear  index 
tn  an  optical  fiber  for  parallel  and  perpendicular  pump  and  probe.  The 
ratio  of  the  two  peaks  is  —3:1. 


The  experiments  agree  with  theory  and  yield  the  correct 
ratio  for  the  two  components  of  n2. 

The  experimental  value  for  n2  can  be  obtained  from  the 
measured  phase  shift.  In  an  experiment,  the  detected  signal 
is  averaged  over  the  pulse  shape  and  the  transverse  area  of 
the  waveguide.  The  measured  peak  of  the  detected  signal  is 

„  SdtSdA  8£pr£ref££„ (2irM)nnLeS 

S~  S  dt  S  dA  4£pr£ref  (13) 

The  detected  signal  is  normalized  to  the  output  of  the  in¬ 
terferometer  when  0bias  =  0  and  Epu  =  0. 

If  the  spatial  field  distribution  in  the  fiber  is  approxi¬ 
mated  by  a  Gaussian,14  and  the  temporal  profile  is  also 
assumed  to  be  Gaussian: 


2ir(2n2||)Z,eff*f  ^llnl/ir 
AtW27T 


(14) 


In  order  to  confirm  that  the  measurement  technique 
functions  as  predicted,  experimental  studies  were  per¬ 
formed  using  an  optical  fiber.  The  optical  fiber  is  an  iso¬ 
tropic  material  with  well-established  nonlinear  optical 
properties.  The  sample  consisted  of  83  cm  of  Coming  non¬ 
polarization  conserving  fiber  with  a  5-jum  core  diameter. 
Special  care  was  required  in  handling  and  mounting  the 
fiber  to  avoid  stress-induced  birefringence  effects. 

In  order  to  demonstrate  the  effect  of  pump  phase,  mea¬ 
surements  were  performed  using  parallel  polarized  pump 
and  probe  pulses,  with  and  without  dithering  the  pump 
phase  <£  ( Fig.  2 ) .  Figure  3  shows  the  experimentally  mea¬ 
sured  phase  shift  for  the  parallel  and  perpendicular  pump 
and  probe  polarization  configurations.  The  two  measure¬ 
ments  were  performed  by  rotating  the  pump  pulse  polar¬ 
ization  90°.  All  other  parameters,  such  as  the  probe  and 
reference  pulse  polarization,  the  phase  bias,  and  the 
intensities  of  all  pulses  were  kept  fixed.  The  pulse  duration 
was  850  fs  and  the  pump  pulse  energy  was  5.5  pJ.  Traces 
were  averaged  over  several  scans  of  the  pump  pulse  delay 
with  a  data  acquisition  time  of  2  min. 

The  ratio  n^/n2L  can  be  determined  from  the  two 
curves.  From  ( 1)  and  (2),  the  ratio  between  the  peaks  of 
the  two  curves  equals  2n2^/rt2l.  The  measured  ratio  is  3.00 
±0.07,  which  implies  n2\/n2L  =  1.50  ±0.04.  This  is  in 
agreement  with  theoretical  predictions  based  on  symmetry 
considerations.  For  an  isotropic  medium,11 

-  (0,0))  =  Xm2(<u:«,  -  (0,(0 ) 

+  Tu2i(<u:^-  -  <o,co)  +  (11) 

In  silica  fibers  the  dominant  contribution  to  the  nonlinear 
index  is  electronic  and  the  three  terms  on  the  right-hand 
side  of  ( 1 1 )  are  approximately  equal.12 13  The  contribution 
from  yields  a  coherent  interaction  between  the  pump 
and  probe  pulses  that  averages  to  zero  when  the  pump 
phase  is  dithered.  Therefore, 

"2|i  _  r(iiii(a>:a>.  ~(Q,(o)  3 

n2i  2x\22\{(0-(0,  —  a>,(o)  2 


where  %  is  the  pulse  energy,  r  the  pulse  FWHM,  and  w  the 
field  radius  (3  jtm  for  the  fiber  used).  Using  the  results  of 
Fig.  3,  the  measured  value  for  «2)j  *s  3.31  ±0.83x  10“  16 
cm2/W.  This  is  in  close  agreement  with  previous  measure¬ 
ments  of  3.2  X  10  ~  16  cm2/W.15  Finally,  it  should  be  noted 
that  the  measurement  technique  is  extremely  sensitive.  The 
signal-to-noise  ratio  in  Fig.  3  is  greater  than  10.  Thus, 
phase  shifts  as  small  as  5  mrad  can  easily  be  detected. 
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We  describe  two-beam  coupling  and  temporal  encoding  experiments  in  barium  titanate.  Volume  gratings  are 
created  in  the  photorefractive  material  by  50-fs  optical  pulses.  Information  in  the  writing  pulses  may  be 
encoded  as  spatially  distributed  volume  gratings  in  the  crystal.  Femtosecond  temporal  waveform  reconstruc¬ 
tion  is  demonstrated. 


The  nonlinear-optical  properties  of  photorefractive 
materials  have  been  extensively  studied  in  wave¬ 
mixing  experiments,  most  of  which  have  been  con¬ 
ducted  using  cw  lasers  as  light  sources.1  Recently 
picosecond  pulsed  lasers  with  low  and  high  repeti¬ 
tion  rates  have  also  been  used  for  the  study  of  photo¬ 
refractive  materials'  *  as  well  as  for  coherence  time 
measurements  of  the  incident  fields.45  The  inter¬ 
action  of  low-energy,  high-repetition-rate  pulsed 
lasers  with  photorefractive  materials  in  the  two- 
beam  coupling  configuration  has  been  studied  theo¬ 
retically,*’  and  the  resulting  coupled  equations  have 
been  solved  in  terms  of  the  coherence  functions  of 
the  incident  fields.  This  theory  describes  multiple- 
pulse  effects  and  the  cumulative  buildup  of  the 
internal  electric  fields  that  are  responsible  for  the 
spatial  modulation  of  the  refractive  index.  It  pre¬ 
dicts  that  photorefractive  properties  are  indepen¬ 
dent  of  the  pulse  duration,  as  long  as  the  pulse 
energy  and  bandwidth  are  not  excessive. 

In  this  Letter,  we  describe  the  nonlinear  interac¬ 
tion  of  femtosecond  pulses  with  barium  titanate 
(BaTiCD  using  a  dispersion-compensated  colliding- 
pulse  mode-locked  (CPM)  laser.  We  demonstrate 
femtosecond  temporal  information  storage  and  ex¬ 
plore  the  factors  that  determine  temporal  resolution. 
Experiments  using  stimulated  photon  echoes"9  have 
demonstrated  the  possibility  of  storing  temporal  in¬ 
formation  in  the  frequency  domain,  with  the  use  of 
inhomogeneously  broadened  atomic  systems.  The 
technique  that  we  present  here  is  more  closely  re¬ 
lated  to  femtosecond  holography.10  In  femtosecond 
holography,  a  spatial  distribution  of  interference 
gratings  results  from  the  different  intersection 
points  of  the  temporal  features  of  the  writing  and 
reference  beams.  This  is  the  approach  taken  by 
Dominic  et  al."  to  perform  autocorrelation  measure¬ 
ments  of  picosecond  pulses  using  photorefractive 
materials.  Conversely,  it  is  also  possible  to  use 
femtosecond  holography  to  convert  spatial  features 
into  temporal  features. 


We  develop  a  simple  geometrical  model  to  describe 
our  experimental  results.  The  model  involves  pa¬ 
rameters  such  as  the  beam  waist  at  the  crystal,  the 
beam  intersection  angle,  and  the  longitudinal  length 
of  the  pulses  (pulse  duration  tp  multiplied  by  the 
speed  of  light).  Similar  models  have  been  used  to 
explain  the  results  of  an  experiment  with  photo¬ 
refractive  materials  in  which  beams  having  partial 
spatial  coherence  were  used12  and  can  also  be  applied 
to  describe  temporal  limitations  of  autocorrelation 
measurements  of  femtosecond  pulses  by  noncollinear 
second-harmonic  generation. 

In  order  to  explore  the  possibility  of  temporal  in¬ 
formation  storage,  we  have  used  the  experimental 
setup  shown  in  Fig.  1.  We  use  a  reference  beam 
and  a  signal  beam.  The  CPM  laser  generates  50-fs 
pulses  at  a  center  wavelength  of  625  nm  and  a  pulse 
repetition  rate  of  100  MHz.  A  Michelson  interfer¬ 
ometer  (M.I.)  is  used  to  generate  a  signal  beam  with 
pairs  of  copropagating  pulses  separated  by  a  vari¬ 
able  time  delay  T.  The  average  powers  used  in  our 
measurements  were  typically  1  mW  in  the  reference 
beam  and  150(x2)  pW  in  the  signal  (double-pulse) 


Fig.  1.  Output  of  a  CPM  laser  is  prechirped  by  a  grating 
pair  and  subsequently  divided  in  two  beams  of  variable 
relative  delay.  A  boxcar  averager,  chopper,  and  autocor¬ 
relator  (AC)  are  used  to  measure  the  duration  of  the 
pulses  after  the  BaTiOi  crystal. 
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Fig.  2.  Longitudinal  length  of  the  pulse  /,,  =  ctp  de¬ 
termines  the  effective  interaction  region  —  defined 
by  Y  =  /sin  f)  —  in  the  limit  of  large  beam  diameter 
non  cos  H  »  Y).  In  the  small-beam-diameter  limit 
icun/cos  H  «  Y ),  the  effective  interaction  region  is  deter¬ 
mined  by  the  intersection  of  the  incident  beams. 

beam.  The  use  of  a  grating  pair  is  essential  in  or¬ 
der  to  precompensate  for  the  material  dispersion. 
Without  dispersion  compensation,  the  pulses  broaden 
to  approximately  200  fs  in  a  single  pass  through  the 
5-mm  BaTiO;t  crystal.  The  beam  diameter  could  be 
varied  by  using  lenses  Ll  and  L2. 

Temporal  information  reconstruction  is  investi¬ 
gated.  A  boxcar  averager  combined  with  a  chopper 
allows  the  pulse,  which  is  scattered  off  the  volume 
grating,  to  be  measured  when  the  copropagating 
signal  beam  is  blocked  by  taking  advantage  of  the 
long  decay  time  of  the  gratings.  The  relatively  long 
dark  period  determined  by  the  chopper  (—1  ms)  also 
ensures  that  free  carrier  gratings  do  not  take  part 
in  the  scattering  process.  It  was  noted  that  the 
buildup  time  and  the  coupling  efficiency  of  the  grat¬ 
ing  under  femtosecond  excitation  is  roughly  the 
same  as  when  a  cw  He-Ne  laser  with  approximately 
the  same  average  power  is  used. 

To  understand  the  role  of  the  geometrical  factors 
that  influence  temporal  information  encoding  and 
reconstruction,  let  us  consider  the  intersection 
region  of  two  beams  in  the  crystal  as  shown  in  Fig.  2. 
If  we  designate  the  wave  vectors  of  the  incident 
beams  k\  and  k  < ,  the  Bragg  condition  is  met  by  light 
scattered  in  directions  k>  and  ku  respectively,  as 
long  as  the  spectral  content  of  the  incident  pulses  is 
small.  The  height  Y  of  the  effective  interaction 
region  in  this  large  beam  diameter  case  is  deter¬ 
mined  by  /,,/sin  0,  where  /„  =  ctr  is  the  longitudinal 
pulse  length.  From  Fig.  2  it  is  clear  that  if  the 
pulse  delay  t  is  changed  bv  a  small  amount,  the  vol¬ 
ume  grating  will  be  formed  at  a  slightly  different  y 
position.  Conversely,  if  the  beam  diameter  is  small 
compared  with  the  longitudinal  length,  the  geomet¬ 
ric  overlap  region  of  the  beams  is  small  and  restricts 
the  region  in  which  the  volume  grating  can  form. 
These  considerations  define  two  limiting  cases 
that  we  explore:  large  beam  diameter  (tun/cos  0  > 
/,/sin  0)  and  small  beam  diameter  (am/cos  0  < 
/,,/sin  B). 

In  order  to  explore  temporal  broadening  effects  in 
the  diffraction  process,  we  first  performed  a  single¬ 
pulse  (one  arm  of  the  Michelson  interferometer 


blocked)  experiment.  The  beam  diameter  at  the 
crystal  was  set  to  ~2  mm.  For  tp  =  50  fs  and 
0  =  3  deg,  u) o /cos  0  >  lp/ sin  0 ,  and  the  mixing  is  in 
the  large-beam-diameter  case  where  the  effective 
interaction  length  is  smaller  than  the  region  of  the 
beam  overlap.  Figure  3(a)  shows  the  autocorrela¬ 
tion  of  the  input  pulses  (after  they  pass  through  the 
crystal),  while  the  autocorrelation  of  the  scattered 
reconstructed  pulse  is  shown  in  Fig.  3(b).  There  is 
a  small  but  noticeable  broadening  in  this  case.  On 
the  other  hand,  when  the  two  lenses  (f  =  15  cm)  Ll 
and  L2  are  introduced,  the  beam  waist  is  reduced 
to  o)0  =  30  pm.  Now  o)0  cos  0  <  lp/ sin  0,  and  the 
mixing  is  in  the  small-beam-diameter  case.  In  this 
case  the  autocorrelation  of  the  scattered  pulses 
shows  almost  no  temporal  broadening  as  seen  in 
Fig.  3(c). 

To  demonstrate  femtosecond  temporal  informa¬ 
tion  storage,  we  use  the  double-pulse  scheme  de- 


Fig.  3.  (a)  Autocorrelation  trace  of  the  prechirped  pulses 
after  a  single  pass  through  the  crystal.  The  autocorrela¬ 
tion  of  the  scattered  pulses  in  the  large-  and  small-beam- 
diameter  limits  is  shown  in  ( b)  and  (c),  respectively.  S.H.. 
second  harmonic. 


Fig.  4  Diagram  showing  that  when  the  temporal  separa¬ 
tion  T  between  the  double  pulses  on  one  of  the  beams  is 
sufficiently  small  relative  to  the  beam  intersection  region, 
temporal  information  encoding  is  possible 
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Fig.  5.  (at  Autocorrelation  of  the  prechirped  double 
pulses  after  a  single  pass  through  the  crystal,  (b)  The 
autocorrelation  of  the  single  pulse  scattered  from  the 
double  gratings. 

scribed  above.  The  situation  that  we  consider  is 
depicted  in  Fig.  4.  When  sufficiently  large-beam 
diameters  are  used,  more  than  one  refractive  index 
grating  can  be  created  in  the  photorefractive  mate¬ 
rial.  In  order  to  encode  information,  the  height  Y 
of  the  effective  interaction  region  must  be  much 
smaller  than  that  of  the  beam-crossing  region  H ,  i.e., 
wu/cos  0  >  cTj sin  0  »  /,,/sin  0,  and  therefore  the 
lenses  LI  and  L2  are  not  used.  For  the  parameters 
of  our  experiment,  the  maximum  pulse  separation 
for  which  two  independent  volume  gratings  can 
be  recorded  in  the  crystal  is  of  the  order  of  T  < 
io n  tan  0/c  ~  1.0  ps. 

In  Fig.  5(a)  we  show  the  autocorrelation  of  the  in¬ 
cident  double  pulses  for  a  temporal  separation  T  of 
450  fs.  In  Fig.  5(b)  the  autocorrelation  trace  of  the 
reconstructed  scattered  pulses  shows  that  the  repro¬ 
duction  of  the  incoming  pulses  is  reasonably  good, 
although  the  wings  of  the  input  pulses  are  more  evi¬ 
dent  in  the  scattered  pulses  and  some  broadening  is 
observed.  From  the  previous  arguments,  it  can  be 
predicted  that  if  the  lenses  are  used  we  would  obtain 
only  a  single  scattered  pulse  because  the  beam  inter¬ 
section  region  is  now  too  small  to  support  two  dis¬ 
tinct  volume  gratings  as  is  in  fact  observed  in  our 
experiments. 

An  interesting  aspect  of  this  double-pulse  experi¬ 
ment  is  that  the  near  field  of  the  scattered  pulses 
(observed  by  blocking  the  double  writing  pulses)  is  a 
spatial  replica  of  the  volume  gratings.  The  near 
field  shows  parallel  bars  produced  by  the  scattering 
of  the  reference  pulse  from  the  stored  gratings. 
The  bars  are  perpendicular  to  the  plane  of  the  inci¬ 
dent  beams,  and  their  separation  changes  in  con¬ 
junction  with  the  temporal  pulse  separation.  These 
bars  are  observed  because  the  intersection  region  of 


the  incoming  beams  is  greater  than  the  crystal 
length  (along  the  propagation  direction). 

In  summary,  we  have  performed  the  first  studies 
of  femtosecond  induced  gratings  in  a  photorefractive 
material.  Temporal  information  encoding  and  re¬ 
construction  have  been  demonstrated.  Pulse  broad¬ 
ening  and  temporal  information  storage  behavior 
are  commensurate  with  a  simple  geometrical  model 
that  uses  spatially  distributed  volume  gratings. 
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We  model  time-resolved  photoluminescence  in  GaAs  using  an  ensemble  Monte  Carlo  method 
coupled  with  a  k*p  calculation  of  the  band  structure.  We  show  that  on  a  picosecond 
scale,  carrier  diffusion  perpendicular  to  the  layer  significantly  reduces  the  density  at  the  surface 
and  consequently  has  a  first-order  effect  on  luminescence  measurements.  To  illustrate 
this  we  compare  the  calculated  luminescence,  with  and  without  diffusion,  to  experimental  data. 


Time-resolved  photoluminescence,  as  recently 
demonstrated,1'5  is  an  excellent  method  for  measuring  car¬ 
rier  relaxation  in  GaAs.  By  measuring  the  time  and  energy 
dependence  of  the  luminescence  with  subpicosecond  reso¬ 
lution,  electron  and  hole  distribution  functions  can  be 
probed  on  the  time  scale  during  which  thermalization  oc¬ 
curs.  By  studying  these  systems,  detailed  information  can 
be  obtained  about  band  structure,  scattering  rates,  and 
nonequilibrium  relaxation  processes. 

Quantitative  analysis  of  time-resolved  luminescence  in 
GaAs  is  difficult,  however,  owing  to  many  different  scat¬ 
tering  mechanisms  and  complex  band  structure.  In  this 
paper,  we  show  that  an  additional  effect,  carrier  diffusion 
perpendicular  to  the  layer,  must  be  included  to  accurately 
interpret  the  luminescence  spectra  on  a  picosecond  scale. 
That  is,  the  density  of  carriers  at  the  surface  decreases  in 
time  due  to  diffusion  into  the  sample.  This  has  a  striking 
effect  on  the  shape  of  time-resolved  luminescence. 

Our  results  show  that  by  omitting  diffusion  in  the  anal¬ 
ysis,  the  carrier  density  is  likely  to  be  overestimated  in 
time-resolved  luminescence  experiments.  This  can  be  seen 
from  the  evolution  of  the  carrier  spatial  profile,  and  from 
the  calculated  distribution  functions  with  and  without  dif¬ 
fusion.  We  show  that  diffusion  has  a  strong  influence  on 
the  luminescence  at  different  energies  and  show  that  it 
must  be  included  to  match  the  experimental  data. 

Our  calculational  approach  uses  an  ensemble  Monte 
Carlo  method6'7  to  model  photoexcited  electron  and  hole 
dynamics.  The  hole  band  structure,  optical  matrix  ele¬ 
ments,  density  of  states,  region  of  photoexcitation,  and 
Bloch  overlap  factors  are  determined  from  a  30-band  k*p 
calculation.8  The  initial  photoexcited  electron-hole  pair 
states  in  the  Monte  Carlo  simulation  are  determined  ac¬ 
cording  to  the  absorption  of  a  monoenergetic  pulse  with  a 
0.5  ps  full  width  at  half  maximum.  The  time-dependent 
distribution  functions  obtained  from  the  Monte  Carlo  sim¬ 
ulation  are  then  used  to  calculate  the  time-dependent  lu¬ 
minescence. 

To  include  diffusion,  the  initial  carrier  spatial  profile 
must  be  determined.  The  Monte  Carlo  method  does  not 
use  a  diffusion  coefficient  because  the  net  redistribution  of 
carriers  depends  directly  on  the  transport  parameters  as 
given  by  the  band  structure  and  scattering  rates.  Thus,  the 


main  ingredient  for  diffusion  in  our  calculations  is  an  initial 
concentration  gradient. 

In  time-resolved  luminescence,  the  photoexcitation 
pulse  is  incident  on  the  surface  of  a  semiconductor  sample. 
Lateral  diffusion  is  unimportant  in  picosecond  mea-  ~e- 
ments  if  the  spot  size  is  chosen  large  enough.  The  spatial 
distribution  of  carriers  perpendicular  to  the  layer  is  given 
by  JVoe_<l0*  for  an  instantaneous  pulse,  where  N0  is  the 
photoexcitation  density,  and  where  a0,  the  linear  absorp¬ 
tion  coefficient,  is  4x  104  cm" 1  for  2.04  eV  photoexcitation 
in  GaAs.9  That  is,  carriers  are  introduced  into  the  Monte 
Carlo  simulation  according  to  the  time  dependence  of  the 
pulse  at  a  k  state  and  position  determined  by  the  absorp¬ 
tion.  For  a  pulse  of  finite  duration,  the  carriers  redistribute 
before  the  pulse  is  completed  resulting  in  a  more  uniform 
profile. 

The  simulated  spatial  profile  of  photoexcited  electrons 
is  shown  in  Fig.  1  for  a  0.5-ps  pulse  with  jV0  =  3X  1018 
cm-3  in  a  4600- A  layer  of  GaAs.  Transparent  confining 
layers  are  assumed.  The  0-ps  curve  in  Fig.  1  is  lowest  be¬ 
cause  at  the  center  of  the  pulse  only  half  of  the  total 
number  of  carriers  are  excited.  After  completion  of  the 
pulse,  the  total  average  density  in  the  layer  is 
N0(  l  _ e~aoL)/a0L,  where  L  is  the  thickness  of  the  layer. 
The  initial  distribution  is  nearly  exponential  because  carri¬ 
ers  have  little  time  to  redistribute  during  a  0.5-ps  pulse. 
Within  several  picoseconds,  however,  the  carriers  diffuse 
away  from  the  surface  and  reduce  the  maximum  density  by 
half. 

As  shown  in  Fig.  1,  the  surface  carrier  density  is 
changing  during  the  rise  time  because  of  diffusion.  It  first 
quickly  increases  with  the  excitation  pulse,  and  then  de¬ 
creases  as  carriers  redistribute.  The  changing  concentra¬ 
tion  complicates  analysis  of  time-resolved  luminescence 
measurements. 

The  expression  for  the  luminescence  is  obtained  by 
summing  the  number  of  photons  over  allowed  optical  tran¬ 
sitions, 

*(«,/)«  i  J  dt'  W(t-f)  1  J  </k  |  //k| 2 

X5(£f(k)-£,.(k)  — <u]/c(k,/')/,.(k,/')  (la) 
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FIG.  1.  Diffusion  of  electrons  into  a  4600-  A  layer  of  GaAs  photoexcited 
by  a  O.S-ps  pulse.  Solid  lines  are  drawn  to  guide  the  eye.  Temporal  center 
of  the  pulse  is  set  at  0  ps.  New  carriers  are  initialized  spatially  according 
to  3 x  1018  e~ cm  \  Because  the  pulse  is  short,  carriers  do  not  redis¬ 
tribute  much  during  the  pulse  and  there  is  still  a  large  gradient  at  1  ps. 
Within  several  picoseconds,  however,  the  electrons  diffuse  away  from  the 
surface  and  reduce  the  maximum  density  significantly. 


cc-  [dr  mt-t')2 1//„(«)|2 

<0  J  U 

Xp„((0  )fc{<0,  f' )/„(«,/'),  (lb) 

where  c  and  v  denote  conduction  and  valence  bands,  H(co) 
is  the  optical  matrix  element,  p(<u)  is  the  joint  density  of 
states,  and  W(t)  is  a  weighting  factor  that  is  determined  by 
the  time  resolution  of  the  experiment.  That  is,  in  the  mea¬ 
surement  the  luminescence  is  sampled  over  a  finite  tempo¬ 
ral  window  given  by  W(t).  A  hyperbolic  secant  squared 
function  with  a  full  width  half  maximum  of  0.5  ps  is  used 
in  these  simulations  for  W(t).  We  assume  collision 
broadening10  and  band-gap  renormalization1 1  are  unimpor¬ 
tant  in  picosecond  measurements  and  neglect  these  effects 
in  our  model.  We  also  assume  that  only  carriers  at  the 
surface  contribute  to  the  luminescence,  since  photons  emit¬ 
ted  deep  in  the  sample  are  reabsorbed  before  reaching  the 
surface.  If  the  sample  thickness  is  much  less  than  the  ab¬ 
sorption  length,  then  all  carriers  can  contribute  to  the  lu¬ 
minescence  and  there  are  no  diffusion  effects.  In  these  ex¬ 
periments,  however,  the  sample  is  more  than  twice  the 
absorption  length  of  0.2  ftm.12 

The  terms  of  most  interest  in  Eq.  1  (b)  are  /c(<u,f)  and 
fv(<o,t),  the  isoenergy  distribution  functions.  These  are  the 
distribution  functions  averaged  over  the  optically  con¬ 
nected  region.  /(&),/)  is  a  function  of  the  luminescence 
energy,  and  is  not  the  same  as  the  energy  distribution  func¬ 
tion,  f(E,t);  if  conduction  and  valence  bands  are  isotropic, 
then  f(E,t)dE  =  f(co,t)d(o.  The  heavy-hole  band  for 
GaAs,  however,  is  highly  anisotropic. 13  The  isoenergy  sur¬ 
face  of  a)  is  determined  by  the  sum  of  electron  and  hole 
energies  for  a  given  k.  A  measure  of  the  time-resolved 
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FIG.  2.  T-vailey  electron  isoenergy  distributions  as  a  function  of  lumi¬ 
nescence  energy  photoexcited  by  a  0.5-ps,  2.04-eV  pulse  (a)  omitting 
diffusion,  and  (b)  including  diffusion.  Without  diffusion  the  maximum 
carrier  density  is  3x  1018  cm  3.  Diffusion  decreases  fe(a>,t )  significantly. 

luminescence,  then,  is  a  measure  of  carrier  dynamics,  albeit 
in  a  nontrivial  manner. 

In  Eq.  (1),  the  recombination  of  electrons  and  heavy 
holes  dominates  the  luminescence.5  Heavy-,  light-,  and 
split-off  hole  dynamics  are  simulated,  but  the  heavy  holes 
control  the  luminescence  because  of  their  overwhelming 
number:  at  1  ps  after  the  center  of  the  excitation  pulse,  the 
Monte  Carlo  simulation  typically  shows  more  than  95%  of 
the  hole  population  is  in  the  heavy-hole  band.  Heavy  holes 
have  very  high  scattering  rates,14  so  convergence  of  the 
distribution  function  occurs  nearly  by  the  end  of  the  laser 
excitation  pulse,  especially  at  high  density  where  hole-hole 
scattering  is  strong.  Electrons  provide  virtually  all  of  the 
structure  on  a  picosecond  scale,  however,  because  they 
cool  much  slower  than  heavy  holes.  Electrons  are  excited 
at  higher  energy  than  the  heavy  holes,  and  inelastic  scat¬ 
tering  rates  are  much  lower.  Electron  relaxation  is  also 
slower  because  intervalley  scattering  prevents  a  direct  de¬ 
scent  to  the  bottom  of  the  T  valley.6 

Since  electrons  provide  most  of  the  structure  of  lumi¬ 
nescence  on  a  picosecond  scale,  in  Fig.  2  we  show  fc(co,t ) 
at  the  surface  in  the  T  valley  for  two  Monte  Carlo  simu¬ 
lations:  (a)  without  diffusion,  and  (b)  with  diffusion.  The 
area  under  /  increases  with  time  because  the  population  in 
the  T  valley  increases  as  carriers  transfer  from  the  satellite 
valleys.  Comparing  Fig.  2(a)  with  Fig.  2(b)  shows  that 
the  occupation  probability  is  less,  and  hence  the  surface 
density  and  the  Fermi  energy  are  less,  when  diffusion  oc¬ 
curs.  As  Fig.  1  shows,  the  surface  density  is  decreasing 
during  the  rise  time  of  the  luminescence.  Diffusion  thus  has 
a  significant  influence  on  the  distribution  functions. 

The  importance  of  diffusion  on  the  luminescence  can 
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FIG.  3.  Calculated  luminescence  in  GaAs  at  1.43,  1.37,  1.66,  and  1.70  eV 
for  the  first  20  ps  after  3xl0‘*  cm-3,  0.5-ps,  2.04-eV  photoexcitation 
compared  with  measurements  by  Block  et  al.  (see  Ref.  1 ).  Experimental 
curves  are  time  shifted,  but  the  time  scale  is  preserved.  Characteristics  of 
the  experimental  data  are  an  increasing  rise  time  for  decreasing  energy. 
The  data  follow  this  pattern  strictly  only  when  diffusion  is  included. 

most  easily  be  seen  by  comparisons  at  different  energies. 
Figure  3  shows  calculated  luminescence  for  3x  101®  cm-3 
photoexcitation  in  GaAs  at  1.45,  1.57,  1.66,  and  1.70  eV 
with  and  without  diffusion,  along  with  time-resolved  lumi¬ 
nescence  measurements  digitized  from  Ref.  1.  Generally, 
rise  time  and  decay  time  of  the  luminescence  decrease  with 
energy.  The  calculations  including  diffusion  agree  better 
with  the  experiment  than  without  diffusion.  This  is  partic¬ 
ularly  evident  at  1.57  eV,  where  without  diffusion  the  lu¬ 
minescence  does  not  decay. 

Understanding  Fig.  3  relies  on  an  explanation  of  the 
decay  of  the  initial  transients.  The  luminescence  does  not 
appreciably  deplete  the  carrier  populations  since  for  these 
experiments  the  recombination  time  is  two  to  three  orders 
of  magnitude  slower  than  the  measurements.  Any  decay  in 
the  luminescence,  therefore,  is  due  to  hot  carriers  relaxing 
to  the  band  edge,  that  is,  passing  through  the  isoenergy 
surface. 

With  this  in  mind.  Fig.  3  can  be  understood  in  terms  of 
Figs.  1  and  2.  Figure  1  shows  diffusion  decreases  the  sur¬ 
face  density  in  the  first  6  ps.  At  1.70  and  1.66  eV  in  Fig.  3 
this  decreasing  density  accelerates  the  decay  when  diffu¬ 


sion  is  included.  At  1.57  eV  without  diffusion,  Fig.  2  shows 
that  fc((o,t)  is  still  increasing  after  20  ps,  thus  causing  the 
luminescence  to  steadily  increase.  With  diffusion  included, 
however,  most  of  the  states  within  a  phonon  energy  below 
1.57  eV  are  unoccupied,  allowing  electrons  to  continue  re¬ 
laxing  and  causing  the  luminescence  to  decay.  At  1.45  eV 
there  is  no  meaningful  difference  between  luminescence 
rise  times,  because  most  states  are  occupied  at  low  energy 
even  with  diffusion. 

We  have  simulated  time-resolved  photoluminescence 
in  GaAs  by  using  a  Monte  Carlo  method  to  model  carrier 
transport  and  by  calculating  the  recombination.  Density  of 
states,  Bloch  overlap  factors,  and  optical  matrix  elements 
were  determined  by  a  k*p  method  to  realistically  include 
hole  band-structure  effects.  Our  results  show  that  including 
carrier  diffusion  in  the  analysis  causes  the  surface  concen¬ 
tration  to  change  on  a  picosecond  scale  and  be  significantly 
lower  than  without  diffusion.  We  compared  our  calculated 
luminescence  with  experimental  data  and  showed  that  dif¬ 
fusion  is  a  critical  consideration. 

We  are  grateful  for  useful  discussions  during  the  de¬ 
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We  present  solutions  to  a  set  of  rate  equations  for  the  electron  dynamics  after  photoexcitation  by  a 
2.0-eV  laser  in  GaAs  and  InP.  Results  obtained,  although  simpler  than  full  Monte  Carlo  solutions, 
closely  follow  the  experimental  data  and  provide  insight  into  intervalley  scattering.  Calculations  show 
that  the  net  return  time  of  electrons  from  the  satellite  L  valleys  into  the  T  valley  is  not  limited  by  the  in¬ 
tervalley  scattering  rate,  but  is  instead  limited  by  the  polar-optic-phonon  scattering  rate  within  the  T 
valley.  This  shows  that  the  time-dependent  mobility  and  luminescence  experiments  depend  on  the  L- 
valley  depopulation  rate,  which  differs  from  the  I  —  T  intervalley  scattering  rate.  Results  further  sug¬ 
gest  that  the  f— -L  scattering  rate  is  faster  than  the  polar-optic-phonon  scattering  rate. 


INTRODUCTION 

Intervalley  scattering  in  compound  semiconductors  is 
responsible  for  the  transferred  electron  effect1  and  is  thus 
important  in  the  operation  of  Gunn  oscillators  and  mi¬ 
crowave  devices.  As  a  result,  there  is  interest  in  accu¬ 
rately  measuring  the  intervalley  scattering  rates.  With 
the  development  of  ultrafast  lasers,  experiments  with  a 
time  resolution  comparable  to  the  intervalley  scattering 
times  have  become  possible.  These  experiments  include 
the  rise  time  of  band-edge  photoluminescence, 3,3  pump- 
infrared-probe  absorption,4  femtosecond  reflectivity,5 
equal  pulse  correlation,6,7  pump-conunuum-probe  ab¬ 
sorption,8  transient  nonlinear  absorpt’  t 0  10  and  cw  ac¬ 
ceptor  luminescence. 11-13 

Several  of  these  experiments  claim  to  measure  the  de¬ 
formation  potential  phonon  scattering  rate  in  GaAs  for 
the  scattering  of  carriers  from  the  satellite  L  valleys  back 
into  the  central  T  valley.  The  values  reported  for  the  de¬ 
formation  potential  constant  DrL,  however,  have  a  large 
variance.13  In  this  paper,  we  show,  through  a  series  of 
rate-equation  models,  that  the  time  measured  in  some  of 
these  experiments  is  not  the  L— -T  scattering  rate,  but 
the  L-valley  depopulation  rate.  This  suggests  that  the  ap¬ 
parent  controversy  is  due  in  part  to  a  misinterpretation  of 
some  of  the  experiments. 

Comparison  of  our  calculations  to  luminescence  and 
mobility  experiments  shows  that  in  GaAs,  the  L-valley 
depopulation  rate  is  most  strongly  influenced  by  the  rate 
of  inelastic  scattering,  chiefly  the  polar-optic-phonon 
(POP)  emission  rate,  within  the  I"  valley.  If  the  inelastic 
scattering  rate  in  the  T  valley  is  small  compared  to  the 
T  — scattering  rate,  then  to  lowest  order,  the  L-valley 
depopulation  rate  does  not  depend  on  the  L  — ►  T  rate. 
For  this  case,  the  bottleneck  for  the  return  of  electrons 
from  the  satellite  L  valleys  to  the  central  T  valley  is  the 
cooling  of  the  electrons  in  the  T  valley.  The  cooling  of 
electrons  allows  a  net  flow  of  electrons  from  the  L  valley 
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into  the  T  valley.  In  addition,  our  studies  show  that  the 
r—L  scattering  rate  is  larger  than  the  POP  scattering 
rate. 

This  paper  is  organized  as  follows.  We  start  with  a 
two-state  rate-equation  model.  While  this  model  is  too 
simple  to  describe  optical  experiments  in  GaAs,  it  illus¬ 
trates  an  important  point,  namely  that  the  two  intervalley 
scattering  rates  T— »L  and  L  —  T  are  not  independent  but 
are  related  through  the  density  of  states  in  each  valley. 
We  then  extend  our  method  to  a  three-state  model.  This 
allows  us  to  describe  actual  experiments  in  GaAs  provid¬ 
ed  that  the  measured  quantities  depend  only  on  the  valley 
the  electrons  occupy  and  not  on  the  details  of  the  elec¬ 
tron  states  within  each  valley.  (An  example  of  such  an 
experiment  is  the  femtosecond  mobility  experiments  of 
Nuss,  Auston,  and  Capasso.5)  Finally,  we  extend  the 
three-state  model  to  a  four-state  model.  With  four  states, 
we  can  account  for  experiments  that  are  sensitive  not 
only  to  which  valley  the  electrons  occupy,  but  also  are 
sensitive  to  which  state  within  the  valley  the  electron  oc¬ 
cupies.  For  example,  in  the  time-resolved  luminescence 
experiments  of  Shah  et  al.,2  only  electrons  at  the  bottom 
of  the  conduction  band  contribute  to  the  luminescence 
signal. 

TWO-STATE  MODEL 

We  start  with  a  two-state  model  shown  in  Fig.  1.  It  is 
the  simplest  possible  model,  but  illustrates  several  key 
points  about  electron  relaxation  in  compound  semicon¬ 
ductors.  One  state  represents  all  electrons  in  the  T  valley 
and  the  other  represents  all  electrons  that  are  in  the  L 
valleys.  yTL  is  the  T  — *L  scattering  rate,  and  yLr  is  the 
L  —  T  scattering  rate. 

This  is  not  a  realistic  model  for  carrier  dynamics  in 
compound  semiconductors  because  it  does  not  allow  for 
states  in  the  T  valley  that  are  not  energetically  able  to 
scatter  into  the  L  valley.  That  is,  relaxation  of  electrons 
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FIG.  1.  The  two-state  model  for  intervalley  scattering  in 
compound  semiconductors.  The  state  T  represents  all  electrons 
in  the  central  T  valley,  while  the  state  L  represents  all  electrons 
in  the  L  valleys.  No  distinction  is  made  between  electrons 
within  a  given  valley.  yrL  is  the  T  to  L  scattering  rate  and  yLr 
is  the  L  to  f  scattering  rate. 


in  the  T  valley  is  ignored.  Nonetheless,  this  model  is  in¬ 
structive  for  the  analysis  of  more  complex  systems  and 
has  features  applicable  to  real  systems. 

The  rate  equations  in  this  model  are  given  in  matrix 
form  by 


"r 

YrL 

YLr 

nr 

di 

nL 

Ytl 

~YLr 

"I- 

Here  nr  and  nL  are  the  densities  of  electrons  in  the  T 
and  L  valleys,  respectively.  A  straightforward  calcula¬ 
tion  shows  that  the  eigenvalues  and  (unnormalized)  eigen¬ 
vectors  are  given  by 


1  0  ,  /*;—  (1'ri  +  Ylt^  * 


YLr 

1 

'  V;  = 

Ytl 

—  1 

(2) 


The  zero  eigenvalue  results  from  the  total  density 
nr  +  being  constant.  The  general  solution  to  the  two- 
state  model  is  obtained  by  taking  a  superposition  of  the 
two  eigensolutions: 


nr(f) 

YrL 

1 

nL[t) 

=  A 

Yll 

+B 

-1 

(3) 


The  constants  A  and  B  are  determined  by  the  initial  con 
ditions: 


Ylt 

1 

A 

f*° 
n  r 

Ytl 

-1 

B 

.0 

nL 

Here  and  are  the  initial  populations  of  carriers  in 
the  two  valleys.  Solving  Eq.  (4)  for  A  and  5,  we  obtain 
the  solution 


nr(f  )  = 


nY  LT 

YrL  ■*" ?'zr 


+ 


nL(t)  — 


nYrL 

YrL+YLr 


nr7'ri  nL.Y  lt 
YrL+YLr 

^h’rL~nLYLr 

YrL+YLr 


-irri+^r*' 

*-  » 

(5) 

g  -l>'n 


The  total  number  of  carriers  n  =  nr(r)  +  ni(t)  =  nr-l-n° 
is  independent  of  time. 

In  equilibrium  (r  — oe  )  the  number  of  carriers  in  each 
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valley  is  found  to  be 


nf  =nr(t-*  oc  )  = 


nYLT 


1*4  = 


nLU-+  oc  )  = 


YrL  +  YLr 
nYn 


(6) 


Yn+YLr 

From  this,  one  obtains  a  “detailed  balance”  relation 

Ylt  _  "r*  _ 

YrL  «r  ~ 


(7) 


R  is  defined  to  be  the  equilibrium  ratio  of  the  populations 
in  the  two  valleys.  If  a  more  detailed  model  based  on  a 
complete  description  of  all  states  were  considered,  then 
intervalley  scattering  events  would  obey  a  detailed  bal¬ 
ance  relation  given  by 

-  (8) 

with  the  transition  probability  per  unit  time  from 

state  k  in  the  T  valley  to  state  k'  in  the  L  valley.14  This 
means  that  the  T— *L  and  L-*T  rates  are  dependent. 
Knowledge  of  one  implies  knowledge  of  the  other. 

Another  point  to  note  is  that  deviations  from  the  equi¬ 
librium  populations  relax  with  a  rate  given  by  the  sum  of 
the  two  scattering  rates  YrL+YLr-  That  is,  if 
8nr(t)  =  nr(t)  —  n fq,  then 


8nl(t)= 


nrYrL~nL.YLr 


Ytl  +  Ylt 

—  6nr(r)  = 


-irn.-’-ru-1' 


nLYLr~nrYrL 


Ytl+Ylt 


(9) 


-'rrz^zr1' 


This  counters  the  naive  assumption  that  a  deviation  from 
the  equilibrium  population  in  the  T  valley  relaxes  with 
yTL  and  a  deviation  in  the  L  valley  relaxes  with  yLr. 
Any  deviation  from  the  equilibrium  ratio  of  carriers  will 
relax  with  the  combined  rate  YrL+YLr-  while  this 
seems  strange  for  electrons  in  the  L  valley  where  the 
£— *T  scattering  rate  is  slower  than  the  T  —*L  rate,  it  is 
important  to  remember  that  fewer  electrons  have  to 
transfer  from  the  L  valley  to  restore  the  equilibrium  ra¬ 
tio.  Equation  (9)  means  that  the  carrier  densities  in  the 
two  valleys  equilibrate  on  a  time  scale  comparable  to  the 
fastest  of  the  two  scattering  times. 


THREE-STATE  MODEL 

The  two-state  model  is  simplistic  since  it  does  not  con¬ 
sider  that  the  T  valley  is  made  up  of  several  states.  In 
particular,  there  is  an  energy  threshold  within  the  T  val¬ 
ley  below  which  electrons  can  no  longer  scatter  into  the 
satellite  L  valleys.  A  more  realistic  model  for  intervalley 
transfer  must  allow  for  this.  Such  a  model  is  given  by  the 
three-state  model  shown  in  Fig.  2.  In  the  three-state 
model,  states  in  the  T  valley  are  separated  above  and 
below  A,  the  energy  threshold  for  transfer  into  the  L  val¬ 
ley.  All  T  valley  electrons  with  enough  energy  to  transfer 
into  the  L  valleys  are  in  the  T*  state;  all  T  valley  elec¬ 
trons  with  energy  less  than  A  are  in  the  T<  state;  and 
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FIG.  2.  The  three-state  model  for  intervalley  scattering.  T> 
is  the  state  that  represents  all  electrons  in  the  T  valley  which  are 
energetically  able  to  transfer  into  the  L  valley  state.  The  T* 
state  represents  all  electrons  in  the  T  valley  that  do  not  have 
enough  energy  to  transfer  to  the  L  valley  state.  Scattering  from 
the  F *  state  to  the  T *  state  occurs  via  inelastic  scattering 
within  the  T  valley  with  a  rate  yA.  This  occurs  mostly  by  POP 
emission. 


electrons  in  the  L  valley  are  in  the  third  state.  Only  elec¬ 
trons  in  state  T  >  can  scatter  into  the  I-valley  state  and 
electrons  in  the  L-valley  state  can  only  scatter  into  the 
T>  state.  Scattering  from  the  T5,  state  into  the  T<  state 
is  given  by  the  scattering  rate  y^  and  occurs  through  in¬ 
elastic  scattering  within  the  T  valley.  This  is  chiefly 
through  POP  emission.13  While  POP  absorption  makes 
scattering  from  below  A  to  above  A  possible,  for  simplici¬ 
ty  it  is  ignored  in  this  model  because  POP  emission 
occurs  more  frequently. 

Note  that  for  simplicity  we  have  not  included  scatter¬ 
ing  into  the  satellite  X  valleys.  For  2.0-eV  photoexcita- 
tion  experiments,  only  a  small  fraction  of  the  photoexcit- 
ed  electrons  are  energetically  able  to  scatter  into  the  X 
valleys  so  that  this  is  a  minor  correction.  For  higher  en¬ 
ergies,  the  X  valleys  must  be  included,  or  the  satellite  val¬ 
ley  state  (and  hence  the  appropriate  scattering  rates) 
modified  to  include  both  L  and  X  valleys.  Inclusion  of  an 
X  valley  should  be  straightforward. 

The  equations  for  the  three-state  model  are  given  by 


«r>(r) 

-IYtl+Y^  Ylt  0 

n  r  >  ( r) 

nLit) 

= 

Yri  ~Ylt  0 

nL(t ) 

«r<(r) 

o 

o 

nr<(r) 

The  eigenvalues  are 

*o=0  , 


(12) 


a.  < 


=  AV0  +  BV~e  ~  +CV0e 


i 


(13) 


with  A,  B,  and  C  determined  by  the  initial  conditions 
0  A.*  k _ 

Q  Tn.*-  Tri*- 

7t.r  +  *~  Yrt."*-*- 

n  ta  y  a 

Solving  for  A,  B,  and  C  and  using  the  relations 
k+k_=yLrYs  , 

<rLr  +  *+)<7i.r  +  *-)=:-rn.7j:r  > 


0 

A 

"r' 

B 

= 

.0 

nL 

C 

.0 

nr< 

(14) 


(15) 


we  can  find  the  final  solution  for  the  populations  of  the 
three  states.  For  the  experiments  we  are  considering, 2,5 
electrons  are  initially  excited  only  into  the  T>  state,  i.e., 
n°  >  ^0,  n£  =0,  n°  <  =0;  then 


nr>(t)=n°r> 
nL[t)  =  n°r> 


7ir  +  *- 

A_f 

k~-k_ 

€ 

€ 

Y  tl 


k^-k. 


i  k-' 
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*-  = 
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(11) 


±ViYrt.+YLr  +  Yj-*YLrr±  j  • 


The  unnormalized  eigenvectors  for  the  matrix  in  the 
three  state  model  are  given  by 


Some  insight  into  the  intervalley  problem  is  obtained 
by  looking  at  the  limiting  form  for  the  eigenvalues  in  Eq. 
(11).  As  mentioned  earlier,  the  detailed  balance  relation 
requires  that  Yvl  and  Ylv  ‘n  the  ratio  of  the  equilibri¬ 
um  populations  R.  Taking  into  account  nonparabolicity 
of  the  T  valley  and  the  fourfold  degeneracy  of  the  L  val¬ 
leys,16  for  2-eV  photoexcitation  in  GaAs,  electrons  are 
photoexcited  0.5  eV  above  the  bottom  of  the  conduction 
band1718  and 
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e<! 
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*0.1  . 


(17) 


This  ratio  is  large  because  the  effective  mass  in  the  L  val¬ 
ley  is  substantially  greater  than  the  effective  mass  in  the 
T  valley.19  We  also  note  that  nonparabolicity  is  an  impor¬ 
tant  consideration  in  determining  this  ratio  since  at  0.5 
eV  nonparabolicity  increases  the  density  of  states  by  60% 
in  the  T  valley.20 

Because  of  the  large  value  of  R,  it  follows  that 
y lt  <<  7rz  anc*  square  root  in  Eq.  (11)  can  be  ex¬ 
panded  to  yield 


A„ 


(7rz  7a  •  > 


~YlvYl 
rn+ra  ‘ 


(18) 


Note  that  for  nonzero  rates  |a*|  » |L_ |- 
There  are  two  interesting  cases  depending  on  how  the 
T— L  intervalley  scattering  rate  compares  with  the  in¬ 
elastic  scattering  rate  in  the  T  valley.  In  the  limit  that 
ya  «7rt>  >•*•.  r —I  intervalley  scattering  is  faster  than 
POP  scattering  within  the  T  valley,  then 


The  three-state  model  can  be  applied  to  a  large  number 
of  experiments  that  are  sensitive  primarily  to  which  val¬ 
ley  the  electron  populates.  An  example  of  such  an  exper¬ 
iment  is  that  of  Nuss,  Auston,  and  Capasso  (NAC).5  In 
the  NAC  experiment,5  the  time-dependent  reflectivity 
was  measured  on  a  femtosecond  scale  to  infer  the  elec¬ 
tron  mobility  as  a  function  of  time.  Since  the  mobility  of 
electrons  in  the  T  valley  is  nearly  independent  of  energy, 
and  the  mobility  of  electrons  in  the  L  valleys  is  negligible 
by  comparison,19  a  measure  of  the  mobility  as  a  function 
of  time  is  a  measure  of  the  density  of  electrons  in  the  T 
valley  (without  regard  to  which  state  within  the  T  valley 
the  electrons  occupy).  In  the  three-state  model,  the  mo¬ 
bility  is  therefore  proportional  to 


M(r)  nr>(r)  +  nr<(r) 

pit-*  oe  ) 


‘r> 


(21) 


To  compare  our  calculations  to  the  experiment,  we 
must  estimate  the  scattering  rates.  Based  on  the  scatter¬ 
ing  rates  of  Schichijo  and  co-workers,23,24  we  estimate  the 
rates  to  be  given  by 


nr  "  r - 

=  — ^r4=-*n-  <i9) 

Yu.  n p 

In  this  case,  A_  does  not  measure  the  L—*T  scattering 
rate,  but  instead  measures  the  product  of  the  inelastic¬ 
scattering  rate  and  the  ratio  of  the  equilibrium  popula¬ 
tions  R.  Even  through  y4  might  be  greater  than  yLr,  the 
net  return  of  electrons  from  the  L  valleys  is  determined 
by  how  fast  electrons  relax  to  lower-energy  levels  within 
the  T  valley,  not  on  how  fast  they  scatter  from  the  L  val¬ 
ley  to  T  valley. 

For  n»}'rt..  *e>  POP  scattering  is  faster  than  the 
T  —*L  rate,  we  obtain 

A-  ~  7'ir  ■  (20) 

This  is,  if  electrons  in  the  T  valley  are  rapidly  scattered 
to  low-energy  states,  then  the  return  time  of  the  L-valley 
electrons  is  limited  by  the  L  — *T  scattering  rate.  Only  in 
this  limit  does  the  L-valley  depopulation  rate  measure  the 
L  — *T  scattering  rate. 

It  is  interesting  to  note  that,  although  the  scattering 
rates  given  by  the  y’s  in  this  model  are  constants,  the 
densities  given  in  Eq.  (16'  are  not  characterized  by  a  sin¬ 
gle  exponential  decay.  Therefore  the  depopulation  rate 
of  a  given  state  [given  by  —  d  \n(n)/dt]  is  time  dependent 
in  spite  of  the  fact  that  the  scattering  rates  are  constant. 
Furthermore,  depending  on  the  values  for  the  scattering 
rates,  the  two  eigenvalues  A+  and  can  differ  by  more 
than  an  order  of  magnitude,  which  leads  to  depopulation 
rates  that  can  change  rapidly  on  a  short  time  scale.  Thus 
the  total  number  of  scattering  events  from  the  T  to  L  val¬ 
ley,  which  is  proportional  to  nr>(r)yri,  changes  rapidly 
with  time  after  the  initial  photoexcitation  in  agreement 
with  earlier  Monte  Carlo  calculations.21,22  This  change, 
however,  reflects  the  time  dependence  of  nr>(f)  and  not 
the  scattering  rate  yri . 


Yri.  =2X  1013  s-1  , 

y  lt 

R  = - =0.1,  (22) 

Yu 

yA  =  5X1012  s_1  . 

In  a  later  section,  we  show  the  sensitivity  of  the  results  to 
these  values.  The  results  of  the  three-state  calculation 
are  plotted  in  Figs.  3-6. 

In  Fig.  3,  the  experimental  mobility  from  the  NAC  ex¬ 
periment5  is  shown  as  the  solid  line.  The  calculated 


time  (p») 


FIG.  3.  The  mobility  vs  time  for  GaAs  as  calculated  from  the 
three-state  model.  The  solid  line  is  the  experimental  work  of 
Nuss,  Auston,  and  Capasso  (Ref.  5),  the  dashed  line  is  the  result 
from  the  three-state  model.  It  is  assumed  that  the  mobility  of 
electrons  in  the  L  valleys  is  close  to  zero  while  the  mobility  of 
the  electrons  in  the  J~  valley  is  independent  of  energy.  The  mo¬ 
bility  is  therefore  proportional  to  the  number  of  electrons  in  the 
T  valley  [«r>(r)  +  nr<(f)].  The  initial  peak  in  the  rate  equa¬ 
tion  solution  near  r  =0  results  from  electrons  initially  photoex- 
cited  in  the  T  valley  that  rapidly  transfer  into  the  L  valley.  This 
occurs  on  an  extremely  fast  time  scale  and  is  not  resolvable  in 
the  experimental  data. 
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FIG.  4.  Sensitivity  of  the  rise  time  of  the  mobility  to  the  in¬ 
elastic  (POP)  scattering  rate  yy  The  solid  line  is  the  experimen¬ 
tal  data  (Ref.  5),  the  dashed  line  the  results  of  the  three-state 
rate-equation  model.  For  the  dash-dotted  line,  the  rate  is  dou¬ 
bled,  whereas  for  the  dotted  line,  the  rate  is  cut  in  half.  As  can 
be  seen,  the  rise  time  of  the  mobility  is  very  sensitive  to  this 
quantity. 


curve  from  the  three-state  model,  Eq.  (21),  is  shown  as 
the  dotted  line.  The  two  are  in  close  agreement  a  short 
time  after  t  =0.  The  peak  in  the  rate  equation  model 
near  f  =0  is  caused  by  electrons  initially  excited  in  the  T 
valley  that  quickly  scatter  to  the  L  valleys.  This  fast  ini¬ 
tial  transient  is  not  resolvable  in  the  experimental  data 
and  indicates  that  the  initial  T— -L  transfer  is  rapid. 

To  test  the  sensitivity  of  the  rate  equation  model,  we 
vary  the  parameters.  In  Fig.  4,  we  show  the  sensitivity  of 
the  rise  time  of  the  mobility  to  the  POP  scattering  rate 
Y±.  The  solid  line  corresponds  to  the  experimental  data 
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FIG.  5.  Sensitivity  of  the  rise  time  of  the  mobility  to  the  in¬ 
tervalley  scattering  time  yrL.  The  solid  line  is  the  experimental 
data  tRef.  S)  and  the  dashed  line  the  calculation.  For  the  dash- 
dotted  line,  the  intervalley  rate  yr!_  is  increased  by  a  factor  of  5 
as  well  as  the  L  to  T  rate  yrr-  thus  keeping  the  ratio  R  con¬ 
stant.  As  can  be  seen  a  fivefold  increase  does  not  significantly 
influence  the  results  indicating  the  general  insensitivity  of  the 
nse  time  to  the  intervailey  scattering.  When  the  rates  are  divid¬ 
ed  by  a  factor  of  5,  the  dotted  line  results.  When  divided  by  a 
factor  of  5,  the  intervailey  scattering  rate  is  slower  than  the 
POP  rate  and  thus  shows  some  discrepancy  with  the  experiment 
at  short  times. 


time  (pi) 


FIG.  6.  Sensitivity  of  the  rise  time  of  the  mobility  to  the 
equilibrium  population  ratio  R.  For  the  dash-dotted  line,  the  L 
to  T  rate  yLr  is  doubled  keeping  yrL  constant  and  thus  dou¬ 
bling  the  equilibrium  population  ratio  R.  For  the  dotted  line, 
the  ratio  R  is  halved.  As  can  be  seen,  the  rise  time  of  the  mobil¬ 
ity  is  very  sensitive  to  this  ratio. 


of  NAC,  and  the  dashed  line  to  the  previous  fit  based  on 
the  rate-equation  model  (cf.  Fig.  3).  For  the  dash-dotted 
line,  the  POP  rate  is  doubled,  whereas  for  the  dotted  line, 
the  POP  rate  is  cut  in  half.  As  can  be  seen,  the  mobility 
and  hence  the  number  of  electrons  in  the  T  valley  de¬ 
pends  strongly  on  the  POP  scattering  rate.  If  the  POP 
scattering  rate  is  low,  ys«yrL,  then  according  to  Eq. 
(19),  it  takes  longer  to  get  a  net  transfer  of  electrons  back 
into  the  (high  mobility)  T  valley. 

In  Fig.  5,  we  show  the  sensitivity  of  the  rise  time  of  the 
mobility  to  the  intervailey  scattering  rates.  To  do  this, 
since  we  have  shown  that  the  L  — >T  and  T  — • L  rates  are 
dependent  [cf.  Eq.  (7)],  we  vary  both  yrL  and  yLT  but 
keep  the  ratio  R  constant.  The  solid  line  corresponds  to 
the  experimental  data5  and  the  dashed  line  the  original 
fit.  For  the  dash-dotted  line,  yTL  and  y  LT  are  increased 
by  a  factor  of  5,  whereas  for  the  dotted  line,  they  are  de¬ 
creased  by  a  factor  of  5.  As  can  be  seen,  the  rise  time  is 
not  strongly  affected  when  the  rate  is  increased,  again 
consistent  with  Eq.  (19),  provided  the  ratio  remains  con¬ 
stant.  When  the  rate  is  divided  by  a  factor  of  5,  however, 
then  one  is  no  longer  in  the  limit  yrL  »yA  and  a  slight 
change  in  the  time-dependent  mobility  is  observed  at 
short  times.  For  slow  intervailey  scattering,  i.e., 
0.2Xyri,  there  is  not  a  rapid  transfer  of  the  initial  elec¬ 
tron  populations  into  the  L  valleys  as  in  the  other  cases. 
If  the  T— *L  intervailey  rate  were  this  slow,  one  should 
see  these  effects  in  the  experimental  data.  The  fact  that 
they  are  not  seen  shows  that  the  T  — *L  rate  cannot  be 
this  slow  and  indicates  that  the  intervailey  rate  yri_  is 
greater  than  the  POP  scattering  rate  y  x 

In  Fig.  6,  we  show  the  sensitivity  of  the  rise  time  of  the 
mobility  to  the  density  of  states  ratio  R.  The  solid  line  is 
the  experimental  data  and  the  dashed  line  is  from  the 
original  calculated  results.  In  all  cases,  we  hold  yrL  con¬ 
stant  and  vary  yLr  to  change  the  ratio  R.  For  the  dash- 
dotted  line,  the  ratio  is  doubled,  whereas  for  the  dotted 
line  it  is  halved.  We  can  see  that  the  data  are  very  sensi- 
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tive  to  this  ratio.  In  fact,  the  curves  are  similar  to  those 
in  Fig.  4,  which  we  should  expect  based  on  Eq.  (19). 

These  results  show  that  the  T— *L  intervalley  rate  is 
faster  than  the  POP  rate.  This  agrees  with  the  transient 
nonlinear  absorption  experiments  of  Rosker,  Wise,  and 
Tang,7  Schoenlein  ex  al and  Becker  ex  al.,910  which 
predict  fast  T— rates.  Also,  even  though  the  L— *T 
rate  is  slower  than  the  POP  rate,  the  return  time  of  car¬ 
riers  from  the  L  valley  is  still  limited  by  the  relaxation  of 
the  electrons  in  the  T  valley. 


FOUR-STATE  MODEL 


The  three-state  model  is  useful  for  determining  the  T - 
and  L-valley  populations  and  is  applicable  for  describing 
experiments  which  depend  only  on  which  valley  the  elec¬ 
trons  occupy.  Other  experiments,  such  as  the  rise  time  of 
the  band-edge  photoluminescence,  are  more  sensitive  to 
the  slate  within  the  T  valley  that  the  electrons  occupy.  In 
particular,  the  photoluminescence  experiments  require 
that  the  electrons  be  at  the  bottom  of  the  conduction 
band  to  contribute.  Therefore  the  three-state  model  is 
not  applicable.  To  account  for  this  additional  structure, 
we  propose  the  four-state  model  shown  schematically  in 
Fig.  7.  In  the  four-state  model,  an  additional  state  is  add¬ 
ed  in  the  T  valley  representing  electrons  at  the  band  edge 
and  is  labeled  T  BE.  The  T  *  state  now  represents  all  the 
electrons  with  energies  ranging  from  one  optic  phonon 
energy  above  the  band  edge  to  energies  just  below  the 
threshold  for  scattering  to  the  L  valley.  For  an  electron 
to  scatter  from  the  top  of  the  T  <  state  to  the  bottom  of 
the  T  valley  into  the  T BE  state,  requires  that  the  electron 
emit  approximately  eight  optical  phonons  (the  first  seven 
emissions  keep  the  electron  within  the  T<  state).  The 
scattering  rate  from  T<  to  I~BE,  denoted  by  yBE,  is  there¬ 
fore  approximately  y±/ 8. 

In  matrix  form,  the  four-state  model  is  given  as 
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FIG.  7.  The  four-state  model.  Here  T>  represents  electrons 
in  the  T  valley  that  are  energetically  able  to  transfer  to  the  L- 
valley  state.  T  *  includes  all  electrons  in  the  T  valley  that  are 
just  below  the  energy  threshold  for  transfer  to  the  L  valley  to 
those  electrons  that  are  just  above  the  band  edge.  I~BE 
represents  the  electrons  in  the  T  valley  at  the  band  edge.  Band- 
edge  photoluminescence  is  sensitive  to  the  number  of  electrons 
in  the  TBE  state.  Scattering  from  T<  to  TBE  involves  the  emis¬ 
sion  of  several  phonons  (approximately  8  in  GaAs)  for  electrons 
to  go  from  just  below  the  threshold  for  intervalley  transfer  to 
the  bottom  of  the  band  edge,  and  the  rate  is  given  by  yBE. 


The  eigenvalues  for  this  matrix  are  easily  found  by  ex¬ 
panding  the  determinant  by  cofactors  of  the  last  column. 
They  are  0,  a+,  and  y  BE. 

One  can  solve  Eq.  (23)  by  noting  that  the  coupled  equa¬ 
tions  for  nr>  and  nL  are  the  same  as  before,  and  that  the 
expression  for  np<  can  be  found  by  using  an  integrating 
factor  and  directly  integrating  the  previous  expressions. 
The  expression  for  n  BE  can  then  be  found  by  integrating 
the  expression  for  nr<. 


nr<(t)  =  y&fo'dt'e  7be"  nnr>(t'), 
"rBE(,)  =  W0'*V<'')  • 


(24) 


From  these  equations,  one  obtains  the  solutions 
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FIG.  8.  Rise  time  for  the  luminescence  in  GaAs  and  InP  as 
calculated  from  the  four-state  model.  The  solid  lines  are  from 
the  experimental  work  of  Shah  et  al.  (Ref.  2).  The  dashed  lines 
are  the  results  of  the  four-state  model.  For  InP,  the  intervalley 
scattering  rates  yrt  and  Fir  are  set  to  zero  keeping  all  other 
rates  the  same.  As  can  be  seen,  the  rate-equation  model  accu¬ 
rately  predicts  the  observed  differences  between  GaAs  and  InP. 
Differences  between  the  experimental  curves  and  rate-equation 
model  at  short  times  originate  from  the  finite  i  mporal  width  of 
the  laser  pulse  not  accounted  for  in  the  four-state  model. 


Results  for  the  four-state  model  are  shown  in  Figs. 
8-12.  The  experimental  luminescence  is  proportional  to 
the  density  of  electrons  at  the  bottom  of  the  band,  i.e., 
the  population  of  T BE.  In  Fig.  8,  we  show  the  experimen¬ 
tal  data  for  the  rise  time  of  the  luminescence  in  GaAs  and 
InP  from  Shah  et  al.2  (solid  lines).  We  also  plot  the  re¬ 
sults  of  the  four-state  model  (dashed  lines).  We  use  a  ygE 
rate  of  6X1011  s-1,  approximately  }  of  F°r  *nP- 
since  the  L  valleys  lie  too  high  in  energy  for  intervalley 
transfer  in  the  Shah  et  al.  experiment,2  we  set  the  inter- 
valley  rates  to  zero  keeping  all  other  rates  the  same.  As 
can  be  seen  from  the  curves,  the  rate-equation  model  ac- 
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FIG.  9.  Sensitivity  of  the  rise  time  of  the  luminescence  for 
GaAs  in  the  four-state  model  to  the  electron-phonon  scattering 
rate  'y^  The  solid  line  is  the  experimental  data  of  Shah  et  al. 
Ref.  2i  and  the  dashed  line  is  the  result  of  the  four-state  model. 
For  the  dash-dotted  line,  the  electron-phonon  scattering  rate  is 
doubled,  while  for  the  dotted  line,  the  rate  is  halved.  As  can  be 
seen,  the  nse  time  of  the  luminescence  is  very  sensitive  to  this 
rate. 


FIG.  10.  Sensitivity  of  the  rise  time  of  the  luminescence  in 
GaAs  to  the  intervalley  scattering  rate.  The  solid  line  is  the  ex¬ 
perimental  curve  of  Shah  et  al.  (Ref.  2),  and  the  dashed  line  the 
result  of  the  four-state  model.  For  the  dash-dotted  and  dotted 
lines,  the  intervalley  scattering  rates  are  changed  so  that  the  ra¬ 
tio  R  remains  constant.  For  the  dash-dotted  line  the  rate  is  in¬ 
creased  by  a  factor  of  5,  while  for  the  dotted  line,  the  rate  is  di¬ 
vided  by  a  factor  of  5.  As  can  be  seen  from  the  figures,  the  rise 
time  of  the  luminescence  is  insensitive  to  the  intervalley  rates 
provided  the  ratio  of  the  equilibrium  populations  is  constant. 


curately  predicts  the  dependence  of  the  rise  time  of  the 
luminescence  for  both  InP  and  GaAs.  Since  the  total 
number  of  electrons  in  the  T  valley  as  a  function  of  time 
does  not  change  between  the  three-state  and  four-state 
models  (only  the  occupation  of  the  different  T  valley 
states  changes),  the  rise  time  of  the  mobility  is  exactly  the 
same  as  before. 

In  Fig.  9,  we  check  the  sensitivity  of  the  luminescence 
rise  time  to  the  POP  scattering  rate.  The  solid  line  is 
from  the  experiment,2  the  dashed  is  from  the  original  fit, 
and  the  dash-dotted  and  dotted  correspond  to  doubling 
and  halving  the  POP  rate,  respectively.  Just  as  for  the 


time  (ps) 

FIG.  11.  Sensitivity  of  the  rise  time  of  the  luminescence  in 
GaAs  to  the  density  of  state  ratio  R.  The  solid  line  is  the  exper¬ 
imental  work  of  Shah  et  al.  iRef.  2)  and  the  dashed  line  the 
solution  to  the  four-state  model.  For  the  dash-dotted  line,  the 
ratio  of  the  populations  is  doubled,  while  for  the  dotted  line,  it 
is  cut  in  half.  The  curves  are  similar  to  those  in  Fig.  9,  showing 
that  the  rise  time  of  the  luminescence  depends  on  the  rate  y  ±  (as 
well  as  y be). 
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mobility,  we  see  that  the  rise  time  of  the  luminescence  is 
very  sensitive  to  this  rate.  The  net  transfer  of  carriers 
from  the  L  to  the  T  valley  depends  strongly  on  the  POP 
scattering  rate  within  the  T  valley. 

In  Fig.  10,  we  show  the  sensitivity  of  the  luminescence 
to  the  intervallev  scattering  rates  keeping  the  ratio  R 
constant.  As  before,  there  is  not  a  strong  dependence  on 
the  magnitude  of  the  rate,  even  when  divided  by  a  factor 
of  5,  provided  the  ratio  is  kept  constant.  We  thus  see 
that  experiments  measuring  the  rise  time  of  the  lumines¬ 
cence  are  insensitive  to  the  intervalley  rates,  and  only  de¬ 
pend  on  the  density  of  states  in  the  two  valleys  and  the 
inelastic  scattering  rate  within  the  T  valley. 

In  Fig.  11,  we  show  the  sensitivity  to  the  ratio  R. 
Again,  the  luminescence  is  sensitive  to  this  ratio  as  pre¬ 
dicted  by  Eq.  (19). 

CONCLUSIONS 

We  have  provided  and  solved  a  set  of  rate  equations  for 
intervalley  scattering  in  compound  semiconductors  such 
as  GaAs  or  InP.  Using  standard  values  for  the  transport 
parameters,  solutions  to  these  rate-equation  models  pre¬ 
dict  the  experimental  trends  in  time-dependent  mobility 
and  luminescence  experiments  quite  well.  Although 
these  rate  equations  are  simpler  than  full  Monte  Carlo 
modeling1718  or  numerical  solutions  to  the  Boltzmann 
equation.15  they  nonetheless  illustrate  several  key  points. 

(i)  The  T— L  and  the  I  — *T  intervalley  scattering  rates 
are  dependent  and  related,  through  detailed  balance,  to 
the  densities  of  states  in  each  valley. 

(ii)  In  time-dependent  mobility  and  time-resolved  pho¬ 
toluminescence  experiments,  one  measures  the  depopula¬ 
tion  rate  of  the  L  valley.  The  depopulation  rate  of  the  L 
valleys  is  not  the  same  as  the  L— *T  intervalley  scattering 
rate  unless  the  inelastic-scattering  rate  in  the  T  valley  is 
significantly  larger  than  the  T  — »L  rate.  From  an  esti¬ 
mate  of  the  transport  parameters  that  are  applicable  to 
GaAs  and  InP,  we  find  excellent  agreement  between  the 


rate  equations  and  experimental  data  for  both  the  time- 
dependent  mobility  and  the  rise  time  of  the  luminescence. 
The  scattering  rates  suggest  that  the  L-valley  depopula¬ 
tion  is  dominated  by  the  energy  relaxation  of  the  elec¬ 
trons  in  the  T  valley.  That  is,  POP  scattering  in  the  T 
valley  acts  as  the  bottleneck  for  the  return  of  the  elec¬ 
trons  from  the  L  valley.  This  return  is  not  limited  by  the 
L  — ►  T  intervalley  scattering  rate.  This  suggests  that  time 
dependent  mobility  and  band-edge  luminescence  experi¬ 
ments  are  not  optimum  for  determining  the  intervalley 
scattering  rates. 

(iii)  The  initial  low  mobility  in  the  NAC  experiment 
can  only  be  explained  by  a  rapid  transfer  of  electrons 
from  their  photoexcited  states  into  the  L  valleys,  suggest¬ 
ing  a  fast  T— *L  rate.  This  is  consistent  with  the  insensi¬ 
tivity  of  the  rise  time  of  the  luminescence  in  GaAs  to  the 
T— *L  rate,  as  well  as  the  transient  nonlinear  absorption 
experiments  of  Rosker,  Wise,  and  Tang,1  Schoenlein 
et  al .,s  and  Becker  et  al.,9,10  which  predict  fast  T —L 
rates.  Thus  this  is  further  evidence  that  the  T  — *L  inter¬ 
valley  scattering  rate  is  faster  than  the  POP  scattering 
rate. 

While  the  rate  equation  models  are  simple,  they 
nonetheless  provide  valuable  insight  into  the  qualitative 
carrier  dynamics  in  compound  semiconductors  and  iso¬ 
late  characteristics  often  lost  in  more  detailed  calcula¬ 
tions. 
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Abstract 

We  report  experimental  and  theoretical  investigations  of  intervalley  scattering  in  Ah 
GaAs.  Femtosecond  absorption  saturation  measurements  are  performed  using  tunable  50 
fs  pulses  and  results  interpreted  using  an  ensemble  Monte  Carlo  simulation  as  well  as  a 
rate  equation  model. 
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Summary 

Intervalley  scattering  in  AlGaAs  has  been  studied  by  a  variety  of  techniques.  Subpi¬ 
cosecond  and  CW  luminescence  spectroscopy!1,2!  and  femtosecond  absorption  saturation!3,4! 
measurements  have  been  used  in  attempts  to  determine  the  intervalley  scattering  rates  and 
deformation  potentials.  While  these  meaurements  were  the  first  to  set  limits  on  the  in¬ 
tervalley  scattering  rates,  they  do  not  directly  measure  the  intervalley  scattering  from  the 
initial  optically  excited  states  and  discrepancies  exist  between  reported  values.  Most  of 
these  investigations  measure  the  return  of  carriers  from  the  satellite  valleys  by  probing  the 
occupancy  of  the  T  valley  at  energies  near  the  conduction  band  edge.  These  measurements 
are  strongly  affected  by  inelastic  scattering  within  the  T  valley!3!. 

We  report  investigations  of  intervalley  scattering  in  AlGaAs  using  tunable  femtosec¬ 
ond  pulses.  By  systematically  varying  the  wavelength  and  spectral  content  of  the  pulses, 
different  scattering  channels  may  be  isolated,  such  as  scattering  to  the  L  satellite  valley. 
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Pump  probe  absorption  saturation  measurement  permits  a  direct  investigation  of  scat¬ 
tering  from  the  initial  optically  prepared  carrier  distributions  and  allows  a  sensitive  and 
direct  measurement  of  intervalley  scattering  processes. 

The  laser  source  in  our  experiments  is  a  dispersion  compensated  colliding  pulse  mode- 
locked  ring  dye  laser  (CPM).  40  fs  pulses  from  the  CPM  are  amplified  in  a  high  repetition 
rate  copper  vapor  laser  pumped  dye  amplifier.  The  amplified  pulses  are  focussed  on  a 
jet  of  ethylene  glycol  to  generate  a  femtosecond  white  light  continuum.  The  continuum 
is  filtered  in  the  frequency  domain  with  a  grating  Fourier  filter!6!  shown  schematically  in 
figure  1.  This  filter  allows  continuous  control  of  both  the  center  wavelength  and  bandwidth 
or  duration  of  the  femtosecond  pulses.  In  addition,  arbitrary  temporal  pulse  trains  may 
be  generated  by  appropriate  masking  in  the  Fourier  plane.  Femtosecond  transient  ab¬ 
sorption  saturation  is  measured  using  pump-probe  techniques  with  differential  detection. 
Experiments  are  performed  using  antireflection  coated  0.2  micron  thick  MBE  samples  of 
AljGaj.jAs  which  are  clad  by  Al0.eGao.4As. 

The  mole  fraction  of  A1  can  be  explicitly  chosen  to  vary  the  band  structure  and  select 
desired  transitions.  For  example,  in  Alo.1Gao.9As,  scattering  to  the  L  valley  is  energetically 
allowed  for  carriers  photoexcited  at  energies  of  ~  1.8  eV  or  greater  (see  Fig.  2).  The  use 
of  frequency  synthesis  techniques  for  these  measurements  represents  a  powerful  approach 
since  the  wavelength,  pulse  duration  and  the  initial  spectral  distribution  of  carriers  can  be 
systematically  controlled.  This  permits  a  comprehensive  investigation  of  carrier  dynamics. 
For  example,  figure  3  shows  pump  probe  measurements  of  transient  absorption  saturation 
behavior  in  Alo.1Gao.9As  at  1.86  eV  and  at  1.72  eV  using  50  fs  pulses.  The  presence  or 
absence  of  the  intervalley  scattering  channel  to  the  L  valley  is  manifest  as  a  change  in  the 
behavior  of  the  initial  femtosecond  transient  in  the  data.  Systematic  measurements  may 
be  performed  as  a  function  of  both  photon  energy  and  wavelength. 

The  determination  of  fundamental  scattering  parameters  directly  from  pump  probe 
data  is  complicated  by  the  fact  that  several  scattering  processes  and  transitions  are  in¬ 
volved.  Thus  we  use  two  theoretical  approaches  to  interpret  our  data.  The  first  technique 
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is  based  on  a  discretization  of  the  Boltzman  transport  equation  and  the  scattering  states 
which  yields  a  simple  rate  equation  model.  This  describes  the  carrier  populations  in  the  L 
valley,  high  in  the  T  valley,  and  at  the  band  edge^.  This  model  represents  a  compact  and 
simple  calculational  approach  which  qualitatively  predicts  the  effects  of  different  carrier- 
carrier  and  intervalley  scattering  rates  on  the  experimental  results.  However,  in  order  to 
obtain  a  more  comprehensive  understanding  of  the  carrier  dynamics  we  have  developed 
a  full  ensemble  Monte  Carlo  simulation  of  the  carrier  dynamics.  40,000  electrons  and 
holes  are  simulated  and  scattering  mechanisms  included  are  electron  -  polar  optic  phonon, 
electron  -  electron  with  time  dependent  screening,  electron  -  intervalley,  hole  -  non-polar 
optical  phonon,  and  hole  -  hole.  The  distribution  functions  generated  by  the  Monte  Carlo 
simulation  are  used  to  predict  differential  transmission  curves  which  may  be  compared  di¬ 
rectly  to  the  experimental  data.  Monte  Carlo  techniques  represent  a  rigorous  approach  for 
studying  carrier  dynamics  since  different  experimental  parameters  can  be  used  as  inputs 
for  the  simulations.  In  addition,  the  effects  of  different  scattering  processes  may  be  studied 
independently  by  isolating  them  in  the  simulation. 

We  gratefully  acknowledge  C.A.  Wang  from  M.I.T.  Lincoln  Laboratory  for  preparing 
the  MBE  samples  used  in  this  investigation.  *LHA’s  permanent  address  is  Dept,  de  Fisica, 
Univ.  Fed.  de  Pernambuco,  Brazil. 
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Figure  Captions 


Figure  1.  Schematic  diagram  of  the  experimental  apparatus.  A  grating  Fourier  filter  is 
used  to  control  wavelength  and  bandwidth. 

Figure  2.  Schematic  bandstructure  of  Alo.1Gao.9As  showing  optical  transitions  above 
and  below  the  L  valley. 

Figure  3.  Transient  absorption  saturation  data  for  pump-probe  energies  above  (1.86  eV) 
and  below  (1.72  eV)  the  L  satellite  valley.  The  photoexcited  carrier  density  is  ~  1018  cm'3. 
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